Realization of a feedback controlled flashing ratchet 
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A flashing ratchet transports diffusive particles using a time-dependent, asymmetric potential. 
Particle speed is predicted to increase when a feedback algorithm based on particle positions is 
used. We have experimentally realized such a feedback ratchet using an optical line trap, and 
observed that use of feedback increases velocity by up to an order of magnitude. We compare two 
different feedback algorithms for small particle numbers, and find good agreement with simulations. 
We also find that existing algorithms can be improved to be more tolerant to feedback delay times. 



Brownian motors 1] induce directed motion of diffu- 
sive particles and are used as models for biological and 
artificial nanomachines. A flashing ratchet [2,131 is a spe- 
cific type of Brownian motor that uses an asymmetric 
ratchet potential V{x), with spatial period L (Fig. [1]) 
that can be turned on (control parameter a — 1) or off 
(a = 0) . If the potential is turned on and off randomly 
or periodically on the time scale of particle diffusion over 
the distance aL (Fig. [ij the mean particle velocity can 
reach up to the order of L/2 per flashing cycle. 

Recent theory work predicts that the performance of 
a flashing ratchet can be dramatically improved by feed- 
back algorithms that control a{t) based on particle dis- 
tributions [4] . Experimental implementation of feedback 
control in ratchets is also predicted to exhibit synchro- 
nization effects [a, |6[ and current reversals [7| , and could 
be used to model paradoxical games [8!| or gating mecha- 
nisms that allow processivity in linear, dimeric molecular 
motors [9|. One specific feedback algorithm maximizes 
the instantaneous velocity (MIV) of overdamped parti- 
cles by setting the control parameter a such that the net 
force on the particles is always zero or positive [J]: 
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Here 8 is the Heaviside function, F{xi) is the force on 
each particle if the ratchet potential were turned on, and 
Xi is the particle position. For small N this scheme re- 
sults in a time-averaged center of mass velocity, <v>, an 
order of magnitude larger than in an optimally operated, 
periodically flashing ratchet [J|. 

Recently, another feedback scheme was proposed that 
predicts a moderate improvement on the performance of 
the MIV scheme for small iV> 1 [10]. This scheme max- 
imizes the net displacement (MND) of the particles that 
is expected when the ratchet potential is turned on: 
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reference point that normally is chosen near the potential 
minimum (2:0 = 0). 

Detailed numerical and analytic studies of feedback 
ratchets have been reported. The MIV scheme was found 
to perform less well than optimized, periodically flashed 
ratchets in the deterministic limit of large A^ [J] , a draw- 
back that can be remedied by introducing a force thresh- 



old [ill, [1^ . Realistic feedback implementations [y, [lO| 
will need to take into account delays [a, IS [7| in the im- 
plementation of feedback, as well as errors [l3| in the 
measurement of particle positions and forces. 

In this Letter, we implement and characterize an ex- 
perimental feedback ratchet for small numbers of diffu- 
sive microbcads in water using optical line traps. We 
implement the MIV scheme and observe the predicted 
increase in <v>, as well as a small additional improve- 
ment using the MND scheme. Our experimental results 
compare well to results from Langevin dynamics simula- 
tions. We find that the MND scheme with feedback delay 
can be further improved by tuning the reference point Xq . 

To create a line trap, a focused laser beam is scanned 
periodically and rapidly along a line such that a trapped 
sphere feels a time averaged optical potential. If the scan 
speed and light intensity are kept constant, the sphere 
diffuses freely along the line, which has been used to 
study the entropic interactions lj| and single file diffu- 
sion of colloids [15] . By modulating the light intensity ap- 




where xq and Xi are measured with respect to the po- 
tential minimum within each potential well, and xo is a 



FIG. 1: Black: a piece- wise linear ratchet potential Vu„{x) 
with period L and asymmetry a = 0.2. Gray: The effective 
potential Vcs{x) (Eq. ([5])) experienced by a sphere of radius 
L/2 trapped in Vun{x) has a reduced asymmetry a^s ~ 0.34, 
and a reduced potential height Vq. 



propriately within each scan, one can create a constant- 
force measuring tool for single molecule experiments [iq , 
or a time-dependent flashing ratchet potential [13| • 

We used an acousto-optic deflector (AOD) to scan a 
1064 nm Nd:YAG laser back and forth at 2 kHz in one 
dimension transverse to the optical axis. To create an op- 
tical trap, we used an oil immersion microscope objective 
(lOOx, 1.4 NA, Leica Model No. 11566014) to focus the 
beam into a microchamber. The 60 ^m deep microcham- 
ber is composed of double sided tape, with a channel cut 
out, sandwiched between a coverslip and a microscope 
slide. Silica spheres, 0.9 /xm diameter, in solution are in- 
jected into the microchamber with a syringe via a drilled 
input hole. The line trap system is used to trap up to 
A'^ = 3 spheres in a 17.4 /im long line containing ten 
ratchet periods [L — 1.74 /im). 

Brightfield images (200x50 pixels) are captured by a 
CCD camera at 130 Hz and are analyzed through Lab- 
VIEW software to find the positions of the spheres, and a 
feedback algorithm (MIV or MND) is then apphed. The 
implementation of this process takes about t = 5 ms per 
frame, much less than the time scale of about 300 ms 
for bead diffusion over a distance aL, which has been 
found to be the upper limit for tolerable feedback delay 
for Af = 1 fl. 

To define a potential shape, a periodic signal is sent to 
the AOD via a function generator. A potential minimum 
(maximum) is created by higher (lower) intensity light. 
However, because of angularly non-uniform AOD output 
and limits in AOD resolution, the resulting output light 
from the AOD will not correspond precisely to the input 
signal. To characterize the flat potential (Fig. [2]), we 
obtained statistics on the free diffusion of particles and 
determined the potential from 
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where k is Boltzmann's constant, T is temperature, N(x) 
is the number of data points at position bin x, and 
Ntotai ~ 10^ is the total number of data points. By 
iteratively measuring the effective potential and making 
adjustments to the input function we achieved a potential 
that was flat within one kT and stable over the course of 
weeks. 

Population sampling in equilibrium (Eq. ([3|)) provides 
little data near the maxima of the ratchet potential when 
Vo is more than a few kT . In those regions we use 



V{x) = --f V {x')dx' 
Jo 



(4) 



where 7 = kT/D is the drag coefficient of the micro- 
sphere {D — 0.44 fiui^/s as measured), and v {x') is 
the ensemble average of the instantaneous velocity of the 
beads at position x' . To characterize the ratchet poten- 
tial we use both Eq. ([3]) and Eq. @i and merge the re- 
sults by matching the slopes near the potential minimum. 




FIG. 2: (a) The measured flat potential (black circles indicate 
experimental error) and the ratchet potential measured using 
Eq. (|3j and Eq. Q (black and gray lines, respectively, also in 
(b)). (b) The measured ratchet potential, averaged over all 
ten periods, yields UaS = 0.35 ± 0.02, and Vo = 10.6 ± 0.3 kT, 
compared to V'cff(a;) (black line) calculated from Eq. ([Sjl, using 
a = 0.2, R = L/2. 



where both methods can be assumed to be reliable. We 
then determine the potential experienced by the beads 
as shown in Fig. [21 where Min with a — 0.2 was used as 
an input to the AOD. 

A finite-size sphere of radius R spatially averages over 
a given one-dimensional ratchet potential V{x) and ex- 
periences the effective potential 
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where S{x') is the cross-sectional area perpendicular to 
the ratchet direction of the bead at position x' . The ef- 
fective potential calculated for our experimental param- 
eters is shown in (Figs. Ill2p . and agrees well with the 
measured potential when V{x') = Vun with a = 0.2 is 
used in Eq. (O and when the result is scaled to match 
the measured Vq (Fig. [2). 

To characterize the performance of our ratchet we must 
allow the possibility of back stepping. For this reason, 
measurements were taken by starting the microsphere at 
the minimum of period 3 (Fig. [H^a)), and the motion of 
the microsphere in response to thermal noise and ratchet 
flashing was recorded until the sphere reached the mini- 
mum of period 1 or 10. Data collection was then halted 
and the microsphere was moved back to the starting po- 
sition. This process was repeated until approximately 
2 X 10^ data points were collected, corresponding to a 
typical total displacement of more than 400 L. Fig. [3] 
shows a sample trajectory obtained by periodically flash- 
ing the ratchet near its experimentally determined opti- 
mal period of 700 ms, resulting in a time-averaged mean 
velocity of <i;>= 0.068 ± 0.01 /xm/s. Using a mirrored 



potential we find <v>= — 0.074±0.01 /xm/s, establishing 
a small negative bias of less than 10 nm/s. 
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FIG. 3: (a) Top black: particle trajectory for feedback- 
controlled ratcheting. Bottom black: a{t) for feedback trajec- 
tory. Gray: particle trajectory and a{t) for periodic ratchet- 
ing, (b) Longer time series of the same trajectories, generated 
by stitching together shorter runs. 

We first present feedback control experiments using a 
potential of strength Vb = 10 kT for the case A^ = 1, 
where the MIV scheme is identical to the MND scheme 
when xo = 0. However, for the MND scheme, it is not 
initially clear whether xq = is the optimal choice, and 
in Fig. |4] we show the <v> obtained using the MND 
feedback algorithm as a function of Xq . We find a max- 
imum <w> of 0.7 /xm/s, an order of magnitude larger 
than the maximum velocity found for periodic flashing, 
in agreement with predictions for A^ = 1 [J, [lO| . 

An initially unexpected finding is that this peak veloc- 
ity occurs for xq « —0.3 /im, 0.15 L to the left of the 
minimum. To better understand the roles that experi- 
mental feedback delay and the precise effective potential 
shape may play in this behavior, we performed simula- 
tions as follows [10|- Particle trajectories were calculated 
by the overdamped Langevin equation 



jx^ = a{t)F{x,{t)) + Ut) 



(6) 



where ^i (t) is a random number chosen from a Gaussian 
white noise distribution that has zero mean and corre- 
lation (C*(i)?i(i')) = 27fcTi5y(5(t - t'), and the external 
force is F{x) = —VV{x). Using for ^(a;) the piecewise 
linear potential Viin(a;) with a = 0.33 and no feedback 
delay, the simulations predict a maximum ratchet veloc- 
ity for Xq — (Fig. [4]), where xo = is always taken at 
the potential minimum. To include a feedback delay r in 
the simulations, a{t) in Eq. ([6]) is replaced by an effective 
delayed value, a{t — r). For t = 5 ms (equal to the ex- 
perimental value), the simulations produce a somewhat 
reduced peak velocity at Xq ~ —0.1 ^m for the piecewise 
linear potential (Fig. \^. When using the same time de- 
lay and V^ff(x) (Eq. ([5])), we find that the peak velocity 
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FIG. 4: (Color online) Average velocity <v> as a function of 
Xo for A'' = 1 and MND. From the top: Simulated data using 
Viin, a = 0.33, r = (solid blue curve); Viin, a — 0.33, r = 
4 ms (dashed green curve); Vea and r = 5 ms (red diamonds); 
and experimental data r = 5 ms (open black circles). 
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FIG. 5: Xq positions (circles) giving maximum <v> from 
curves such as the inset of Fig. [6] for various r. For small r 
these points fit well to a diffusive curve {x — \/2Dt, black 
line). Gray line: x = —Fr/'y. 



is shifted even further to negative xo, to approximately 
the same position as in the experiment. 

To understand why finite feedback delay can introduce 
a shift of the maximum <v> to negative xq values, con- 
sider a situation where the potential is on and a particle 
is drifting down the long slope of the ratchet. To achieve 
optimal <v>, the potential should turn off at the in- 
stant when the particle reaches the potential minimum 
y . For finite r it is therefore advantageous to anticipate 
the arrival of the particle at the potential minimum, and 
to trigger the potential to be turned off at a position that 
is passed by the particle on average at a time r earlier. 
To test this explanation, we show in Fig. [5] the Xq that 
produces maximal <v> as a function of t from simula- 
tions using V\in{x) with a — 0.33. Also shown are the 
distance xq — —Ft/^ a particle on the long slope would 
drift during r, and the average distance xq — \J2Dt a 
particle would diffuse during r. For short r, where trans- 
port is mainly diffusive, the optimal choice of xo is indeed 
closely approximated by the diffusive prediction. 

An important implication of Fig. 5] is that existing 
feedback algorithms can be improved to better tolerate 
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ing devices akin to Maxwell's demon, or to study particle- 
particle interactions in ratchets. 
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FIG. 6: (Color online) Simulated data using Vun and a = 
0.33. Black diamonds are from the MIV and MND scheme 
for xo = 0. Open red circles are the maximum <v> obtained 
using MND and by tuning xq for each value of r. Inset: <v> 
as a function of xo for r = to 0.005 L /D (top to bottom) 



an implementation delay r, by introducing a spatial trig- 
ger Xo < that anticipates the arrival of particles at the 
potential minimum. Fig. [S] illustrates this for the case 
A^ = 1. The lower curve shows the data reported in the 
seminal papers p, |6[ for the MIV scheme (equivalent to 
the MND scheme for xq = 0). The upper curve shows 
the optimal <v> found by simulating the MND scheme 
and varying xo (see inset), demonstrating a substantial 
improvement for finite r. In simulations not shown here 
we find that for iV > 1 the optimal MND scheme requires 
a finite xo even when r = 0. A future systematic evalu- 
ation of using spatial triggers in feedback algorithms for 
varying A'' and t is likely to lead to further improvements 
over existing algorithms. 

For small A^ > 1 it has been predicted that the MND 
scheme can yield a moderately higher < v > than the 
MIV scheme [10|. To test this prediction, we performed 
experiments for A^ = 2,3. To avoid any effects from 
particle-particle interactions, data points recorded when 
more than one bead was within any given ratchet period 
were excluded. For A^ > 3 this happens so frequently that 
sufficient data cannot be collected. In Fig. [7ja) we show 
experimental results for the ratio of average velocities 
obtained using the two schemes. In each case, xo was 
varied, and the maximum <v> was used in Fig [TJ^a). In 
agreement with simulations (Fig [7jb)) that used Vim (a;) 
with a = 0.33, t = and optimal xq for MND, the MND 
scheme performs better than the MIV scheme for A^ = 3 
and Vo > 30fcr. A difference not shown here is that 
in the simulation <v> increases with Vq in all cases, 
whereas in the experiment we see an unexplained overall 
drop between Vq = 30/cT and AOkT. This may be due 
distortions in V{x) at high Vq. 

We have demonstrated diffusive particle motion in cus- 
tomized, time dependent potentials with fast feedback. 
In the future this system can also be used to realize sort- 
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FIG. 7: Ratio of <v> obtained for MND and MIV methods 
for various N and Vo = lOfcT, 30fcr, 40fcT from bottom to 
top. (a) Experimental data for optimal xq. (b) Simulated 
data using Viin, a = 0.33, t — and optimal xq for MND. 



